EUROPEAN LABORATORY FOR PARTICLE PHYSICS 



CERN-EP/99-163 
February 7, 2008 



Bose-Einstein correlations in K ± K ± pairs 
from Z° decays into two hadronic jets 

The OPAL collaboration 

Abstract 

Bose-Einstein correlations in pairs of charged kaons produced in a sample 
of 3.9 million hadronic Z° decays have been measured with the OPAL 
experiment at LEP. Charged kaons were identified in the central tracking 
detector using their specific energy loss in the drift chamber gas. The 
correlation function was studied in two-jet events using a double ratio, 
formed by the number of like-sign pairs normalised by a reference sample 
in the data, divided by the same ratio in a Monte Carlo simulation. The 
enhancement at small values of the four-momentum difference of the pair 
was parametrised using a Gaussian function. The parameters of the Bose- 
Einstein correlations were measured to be A = 0.82 ± 0.22 ^ q]^ for the 
strength and Rq = 0.56 ± 0.08 1 ggg fm for the kaon source radius, where 
the first errors are statistical and the second systematic. Corrections for 
final-state interactions are discussed. 
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1 Introduction 



Intensity interferometry was applied in 1953 by Hanbury-Brown and Twiss |T[] in 
radio astronomy in order to estimate the spatial extension of stars (HBT effect). 
In particle reactions which lead to multi-hadronic final states the HBT effect 
manifests itself as a constructive interference between two identical bosons - 
the so-called Bose-Einstein (BE) correlation, which is now well known and was 
first observed by Goldhaber et al. |2| in pp annihilations. There is interest in 
the quantum mechanical aspects of the BE correlations, but they are also used 
to estimate the dimensions of the source of the identical bosons. BE correlation 
studies of pion pairs have been carried out in a large variety of particle interactions 
and over a wide range of energies 0-§f . Recently, pion BE correlations have been 
investigated in connection with the W mass measurement in e + e~ — > W + W~ 
reactions at LEP at centre-of-mass energies above 161 GeV ||]7]]. 

Compared to the abundant information now available on BE correlations in 
pion pairs, knowledge of the correlations in identical strange boson pairs is scarce 
and, until recently, was mainly limited to the KgK° system. K°Kg pairs may 
exhibit a BE correlation enhancement near threshold even if the origin is not a 
K°K° or K° K° system but a K°K° boson-antiboson system. A K°Kg low-mass 
enhancement has recently been observed in hadronic Z° decays ft-P]. However, 



the interpretation of this enhancement as a pure BE correlation effect is compli- 
cated by the possible presence of the f (980) decay into K°K°. Recently it has 
been pointed out that the information from BE correlation studies of the K ± K ± 
system, which cannot result from the fo(980) decay, can serve as an effective tool 
in setting a limit on the resonant fraction of the K°Kg BE enhancement |12 . 



The study of K ± K ± BE correlations also has a bearing on the relation between 
the dimension of the emitting source and the mass of the emitted particles; from 
recent measurements it has been pointed out that the source dimension seems 
to decrease as the mass increases |Tj|. Several models have been proposed to 
account for this behaviour |13| , [14]] . 

This paper presents a study of K ± K ± BE correlations using the high statistics 
sample of Z° hadronic events recorded by the OPAL detector at LEP. The paper 
is organised as follows. In Section || the methodology used for measuring the 
BE correlations is presented. The event and track selections are described in 
Section [3|. In Section f| the analysis of the data is presented, and in Section |5| the 
systematic effects are studied. Finally, the conclusions are drawn in Section ^|. 

2 Analysis method 

The BE correlation function for two identical bosons is defined as: 

p{Pl)P{P2) 



4 



where p(pi,P2) = (l/°")(d 8cr /d 4 Pid 4 P2) is the two-particle phase space density 
subject to BE symmetry, and p(p) is the corresponding single particle quantity 
for a particle with four-momentum p. The correlation function can be studied as a 
function of the four-momentum difference of the pair, Q, where Q 2 = — (pi — P2) 2 
= M 2 — 4m^ oson , and M is the invariant mass of the pair of bosons each of mass 

^■boson • 

From the study of the correlation in pairs one can determine the geometrical 
and dynamical properties of the emitting source. For a static sphere of emitters 
with Gaussian density, the correlation function is parametrised with the Gold- 
haber function || as: 

C(Q) = 1 + Ae-^ 2 , (2) 

where A, the strength of the correlation, is for a completely coherent source and 
1 for a completely incoherent one. The parameter R in GeV -1 is related to the 
radius of the source, Ro, through the relation R = RTic. 

The two-particle phase space density p(pi,P2) is obtained from a sample of 
pairs of identical bosons. In this analysis this sample, called the like-sign sample, 
is formed by pairs of kaons with the same charge. The denominator of Equation |] 
p(pi)p(p 2 ) is in practice replaced by a reference distribution Po(pi,P2), which 
resembles p(pi,P2) in all aspects except in the BE symmetry. A perfect reference 
sample should have the following properties: absence of BE correlations, and 
presence of the same correlations as in the the like-sign sample, arising from 
energy-momentum and charge conservation, the topology and global properties 
of the events, and resonance or long-lived particle decays. 

The principal difficulty for measurements of BE correlations is in the defi- 
nition of the reference sample from which po{pi,P2) is obtained. When corre- 
lations amongst like-sign boson pairs are measured, the obvious reference sam- 
ple is provided by unlike-sign boson pairs. Unfortunately, the Q distribution of 
the unlike-sign pairs includes prominent peaks due to neutral meson resonances 
(e.g. cf) — > K + K~). An alternative reference sample can be derived from Monte 
Carlo simulations in which BE correlations are not included. This relies on a 
correct simulation of the physics in the complete absence of any BE effect, and 
a correct modeling of detector effects. A third type of reference sample can be 
obtained using the methods of event- or hemisphere-mixing, where particles from 
different events or hemispheres are combined. The last, least model-dependent 
method was used in the present work. 

Each event was divided into two hemispheres separated by the plane perpen- 
dicular to the thrust axis and containing the interaction point. A charged kaon 
track in one hemisphere was combined with a kaon track of the same charge found 
in the opposite hemisphere after reflecting the momentum of this second track 
through the origin. To ensure that the like-sign and reference samples were inde- 
pendent, the two kaon tracks forming a pair in the like-sign sample were required 
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to be in the same hemisphere. If all the events have two back-to-back jets, the 
reference sample of hemisphere-mixed pairs will be similar to the like-sign pair 
sample apart from the lack of BE correlations. Monte Carlo studies have shown 
that the hemisphere-mixing technique only works in symmetric topologies such 
as back-to-back jets. Therefore, two-jet events were selected by requiring a high 
value of the thrust. The event- and hemisphere- mixing techniques have already 
been used in other experiments |]15|,|16||. 



3 Event and track selection 

A detailed description of the OPAL detector can be found in reference [[H]]. This 
analysis is primarily based on information from the central tracking detectors, 
consisting of a silicon microvertex detector, a vertex drift chamber, a jet chamber 
and z-chambers[J all of which lie within an axial magnetic field of 0.435 T. The jet 
chamber, which has an outer radius of 185 cm, provides up to 159 measurements 
of space points per track with a resolution in the r-0 plane of about 135 fim and 



a transverse momentum^ resolution of <J Pt /pt = y (0.02) 2 + (0.0015pr) 2 , with p T 
in GeV/c. Particle identification in the jet chamber is possible using the measure- 
ment of the specific energy loss dE / dx |1| with a resolution of approximately of 



3% for high-momentum tracks in hadronic decays Since the identification 

of charged kaons using dE/dx is crucial to this analysis, the calibration of the 
energy loss over the many years of data taking was checked and improved when 
necessary. Control samples of particles identified by techniques other than the 
energy loss were used to remove year-to-year variations in the measured dE/dx 
and to resolve discrepancies between the measured dE/dx in the data and the 
theoretical dE/dx in the Monte Carlo. 

This analysis used a sample of about 3.9 million hadronic Z° decays recorded 
at LEP between the years 1992 and 1995. A sample of 6.75 million Monte Carlo 
hadronic events generated with JETSET 7.4 J2(| and tuned to reproduce the 



global features of the events |21] was also used. The generated events were pro- 
cessed through a detailed simulation of the experiment ]22| and subjected to the 
same event and track selection as the data. A detailed description of the selec- 



tion of hadronic events is given in |23fl. Events with two clear back-to-back jets, 



necessary for the proper functioning of the hemisphere-mixing technique, were 
selected by requiring that the thrust value was larger than 0.95. About 30% of 
the events passed the thrust cut. 

Charged tracks were required to have a minimum transverse momentum of 
150 MeV/c, a maximum reconstructed momentum of 60 GeV/c, a distance of 

1 A right handed coordinate system is used, with positive z along the e~ beam direction and 
x pointing towards the centre of the LEP ring. The polar and azimuthal angles are denoted by 
and (j>, and the origin is taken to be the centre of the detector. 

2 The projection onto the plane perpendicular to the beam axis. 
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closest approach to the interaction point less than 0.5 cm in the plane orthogonal 
to the beam direction and the corresponding distance along the beam direction of 
less than 40 cm. The first measured point had to be within a radius of 75 cm from 
the interaction vertex. The cuts on the transverse and longitudinal distances of 
closest approach help to remove particles from long-lived decays. About 50% of 
the reconstructed tracks passed these selections. 

Kaons were identified using the dE/dx measurement. Only tracks with at 
least 20 hits available for the measurement of dE/dx and with a polar angle 
satisfying |cos#| < 0.9 were considered. For each track, a \ 2 probability was 
formed for each stable particle hypothesis: e, fi, n, K and p. A track was identified 
as a kaon candidate if it had a probability of at least 10% of being a kaon and if 
the probability of being a kaon was larger than the probability of being any other 
of the above particle species. In addition, in order to reject pions, each track 
was required to have a pion probability less than 5%. Electrons from photon 
conversions were rejected using a neural network algorithm as described in [ [24] 



According to the Monte Carlo, approximately 34% of the true kaons passed these 
requirements and the estimated kaon purity of the track sample was about 72% 
on average, with variations between 50% and 97% depending on the momentum 
of the track. The lowest purity corresponded to the momentum range between 



0.9 and 1.5 GeV/c where the pion and kaon bands overlap in dE/dx ||19|| . In the 



selected kaon track sample the fraction of pions was estimated as 17%, that of 
protons as 11%, and the contribution of muons and electrons was negligible. 



4 Data analysis 

The Q distributions of the like-sign kaon candidate pairs and the hemisphere- 
mixed kaon candidate pairs were determined using tracks that passed the selec- 
tions described in Section |[ In the data, 76063 like-sign and 98558 hemisphere- 
mixed kaon candidate pairs with Q < 2.0 GeV were selected. In the Monte Carlo, 
the corresponding numbers of kaon candidate pairs were 109601 and 136785 re- 
spectively. 

Since the BE correlation manifests itself only for identical particles, it was 
necessary to correct for impurities. The Monte Carlo predicted that the main 
contamination in the sample of like-sign kaon candidate pairs was from K7r and 
Kp pairs. In this sample, and for pairs with Q < 2.0 GeV, the estimated fraction 
of KK pairs was about 48%, the K-zr fraction about 27% and the Kp fraction about 
13%. The contamination from pairs susceptible to BE correlations or Fermi-Dirac 
anti-correlations (i.e. im and pp) was negligible, of the order of 3% for nn and 
less than 1% for pp pairs. The fraction of KK pairs was constant over the whole 
Q range. The sample of hemisphere-mixed kaon candidate pairs had an almost 
identical composition. 

The Monte Carlo non-KK Q distribution was subtracted from the data distri- 
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bution using the fraction given by the simulation and normalised to the number 
of data pairs. Both like-sign and hemisphere-mixed pair distributions were cor- 
rected using this technique. Figure [I] shows these corrected Q distributions for 
the data and the Monte Carlo in which BE correlations were not simulated. The 
hemisphere-mixed Q distributions were normalised to the like-sign Q distribu- 
tions in the region 0.6 < Q < 2.0 GeV where no BE correlations are expected. 
Both data and Monte Carlo distributions show a similar behaviour at high val- 
ues of Q: the hemisphere-mixed distribution is above the like-sign one in the 
region 0.7 < Q < 1.2 GeV; there is a cross-over of both distributions at Q about 
1.2 GeV; and the like-sign distribution is above the hemisphere- mixed one at 
Q > 1.2 GeV. 

Figure |^ shows the ratio N ++ {Q) / N mix {Q) in both the data and the Monte 
Carlo, where N ++ (Q) is the number of like-sign pairs and N mix (Q) is the number 
of hemisphere-mixed pairs as functions of Q. The first two bins of the distri- 
bution shown in the figure were combined due to the small statistics. The data 
distribution shows a clear enhancement in the region Q < 0.3 GeV. There is also 
a rise of the correlation at high values of Q normally attributed to long-range 
correlations. Indeed, the slope of the correlations at high Q is well reproduced 
by the Monte Carlo. The ratio N ++ (Q) / N m i X (Q) in the Monte Carlo deviates 
slightly from unity at high Q and falls slowly with decreasing Q, probably due to 
features of string fragmentation and local conservation of charge and strangeness. 
These effects can be taken into account if both like-sign and hemisphere-mixed 
data distributions are divided by the corresponding Monte Carlo distributions. 

The correlation function was therefore defined as the double-ratio: 

r Nga^ N MO( Q) 

mix[Q) " N^{Q)I N%g{Q)> [6) 
and was parametrised using a modified version of Equation ^ 

C(Q) = N(1 + Ae-W* 2 ) (1 + 5Q + eQ 2 ), (4) 

where N is a normalisation factor, and the empirical term (1+5 Q+eQ 2 ) accounts 
for the behaviour of the correlation function at high Q values due to any remaining 
long-range correlations. Figure |3] shows the correlation C m i X (Q) with the result 
of the fit. The fitted parameters and the correlation coefficients between them 
are given in Table |l|. The fit has a % 2 = 43 for 35 degrees of freedom. 



5 Systematic effects 

Systematic effects arising from the event and track selections, the modeling of 
dE/dx in the Monte Carlo, the parametrisation of the correlation function and 
the choice of the reference sample are considered. In each case, the result of the 
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Correlation coefficients 




fitted value 


N X R 5 e 


N 
X 

Ro (fm) 
5 (GeV- 1 ) 
e (GeV" 2 ) 


0.97 ± 0.11 
0.82 ± 0.22 
0.56 ± 0.08 
-0.07 ± 0.16 
0.06 ± 0.06 


-0.59 +0.52 -0.96 +0.85 
+0.17 +0.58 -0.52 
-0.42 +0.29 
-0.96 



Table 1: Fitted parameters and correlation coefficients obtained in the Gaussian 
parametrisation of the correlation function C m i x {Q). The uncertainties on the 
parameters are statistical only. 



fit to the correlation function is given in Table 0, with row (a) giving the result 
of the basic fit discussed in the previous section. 

The overall systematic uncertainties in the parameters A and Ro were cal- 
culated as the largest single deviations between the parameters of the fits from 
rows (b) to (1), and the parameters of the basic fit in row (a). The final values of 
the parameters are 

A = 0.82 ± 0.22 t ° \l 
R = 0.56 ± 0.08 t aoe fm. 

The mismodeling of the kaon momentum spectrum, residual BE correlations in 
the reference sample, the origins of the kaons and final-state interaction correc- 
tions are also discussed in this section. 



5.1 Event and track selection 

• As discussed in Section ||], the hemisphere-mixing technique only works 
when two-jet events are selected. The analysis was repeated for events 
selected using a cone jet finding algorithm p5f instead of the cut in thrust 



(row (b)). The value of the thrust cut was also changed from 0.95 to 0.93 
and the analysis repeated (row (c)). 

To obtain a purer sample of kaons, tracks with momenta in the pion-kaon 
dE/dx overlap region, 0.9 < p < 1.5 GeV/c, were rejected (row (d)). 

The minimum number of dE/dx hits required for each track was increased 
from 20 to 40 (row (e)). 



5.2 Parametrisation of dE/dx 

Since the analysis relies to a large extent on the separation of pions from kaons, it 
is especially important to understand the dE/dx measurements of the copiously- 
produced pions. A sample of pions was identified in K° — > tt + tt~ decays || 
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and was used to estimate the mismodeling of the normalised ionisation energy 
loss f24|. The normalised dE/dx is defined as ^[gy^ = (dE/dx — 

(dE /dx) ) I \o{dE /dx)) . Here, dE/dx is the measured value, while (dE/dx) 
and (a(dE/dx)) are the expected value and the expected error assuming the 
track to be a pion. This analysis showed that the mean value and the width of 
the normalised dE/dx distribution were both known to within ±10% of o. By 
changing in the Monte Carlo the normalised dE / dx of both kaons and pions by 
their known uncertainties, the fraction of non-KK pairs was found to vary by up 
to ±10%. The BE analysis was consequently repeated with the assumed impurity 
value set to 57% and 47% (rows (f) and (f)). 



5.3 Fit of the correlation function 

• The binning of the Q distributions was modified from 50 MeV to 20 MeV 

(row (g)). 

• The Q distribution normalisation range was changed from 0.6 <Q< 2.0 GeV 
to 0.8 < Q < 2.0 GeV (row (h)). 

• The fit was repeated with the first bin Q < 0.05 GeV, excluded. This was 
done to test the effect of potential problems at low Q because of the limited 
resolution p6f (row (i)). 

• If the use of the double-ratio removes all correlations other than BE, such as 
long-range correlations at high Q, then the empirical term (1 + 8Q + eQ 2 ) 
of Equation (| would not be necessary. C m i X {Q) was parametrised using the 
simplified function: 

C(Q)=N(l + \e- {RQ)2 ). (5) 

The parameters of the fit are given in row (j). The fit has \ 2 = 49 for 37 
degrees of freedom . This fit was also repeated with the range limited to 
Q < 1.5 GeV to reduce possible long-range correlations (row (k)). 



5.4 The reference sample 

The BE correlation was also measured using Monte Carlo like-sign pairs as the 
reference sample in the sub-sample of two-jet events (row (1)). The correlation 
function in this case was defined as: 

N data (Q) 

Cmc{Q) = J^McJ^- (6) 

Figure f| shows the correlation function Cmc(Q) together with the results of a fit 
using Equation f|. The values of the parameters are: iV = 0.88 ± 0.05; A = 0.92 
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± 0.17; R = 0.59 ± 0.06 fm; 5 = 0.04 ± 0.08 GeV" 1 ; e = 0.05 ± 0.03 GeV~ 2 , 
where the errors are statistical only. The fit has x 2 = 42 for 35 degrees of 
freedom. 

Although there are known imperfections in the simulation, particularly at low 
momenta in the kaon momentum spectrum (see section |5.5| ), these results were 
taken as an indication of systematic differences between the hemisphere-mixing 
method and the simple use of a Monte Carlo reference sample. 



fit variation 


A 


#o (fm) 


xVDoF 


(a) basic fit 


0.82 ± 0.22 


0.56 ± 0.08 


43/35 


(b) cone jet finding 


0.90 ± 0.21 


0.58 ± 0.13 


45/35 


(c) thrust > 0.93 


0.88 ± 0.19 


0.53 ± 0.06 


42/35 


(d) cut on p 


0.89 ± 0.23 


0.58 ± 0.08 


37/35 


(e) minimum N^E/dx — 40 


0.79 ± 0.27 


0.62 ± 0.13 


43/35 


(f) dE/dx +10% 


0.99 ± 0.26 


0.55 ± 0.07 


44/35 


(F) dE/dx -10% 


0.70 ± 0.20 


0.58 ± 0.11 


42/35 


(g) Q binning of 20 MeV 


0.77 ± 0.22 


0.50 ± 0.06 


80/95 


(h) normalisation 0.8 < Q < 2.0 GeV 


0.82 ± 0.21 


0.56 ± 0.08 


43/35 


(i) Q lower limit 0.05 GeV 


0.82 ± 0.23 


0.57 ± 0.09 


43/34 


(j) fit using Equation |5] 


0.85 ± 0.21 


0.64 ± 0.08 


49/37 


(k) Q upper limit 1.5 GeV 


0.81 ± 0.21 


0.59 ± 0.08 


38/28 


(1) Monte Carlo as reference sample 


0.92 ± 0.17 


0.59 ± 0.06 


42/35 


Overall systematic error 


+ 0.17 
- 0.12 


+ 0.08 
- 0.06 





Table 2: Results of various fits of the BE correlation function; the quoted 
errors associated with A and Ro are statistical only. The results in row (a) 
correspond to the basic fit of Section f| while the other lines show the results 
obtained when modifying certain criteria as explained in the text. The overall 
systematic uncertainties in the parameters A and Rq were calculated as the largest 
single deviations between the parameters of the fits from rows (b) to (1), and the 
parameters of the basic fit in row (a). The final systematic uncertainties in A and 
Ro are given in the last row. 



5.5 Mismodeling of the kaon momentum spectrum 

If the simulation were perfect, one would expect to get the same results by mea- 
suring the correlation function with C mix (Q) and with Cmc(Q)- The double-ratio 

r N d + t(Q) jgg(g) _ N^{Q) JNg(Q) 

N mix(Q) N mix{Q) ' N ++(Q) N mix(Q) ^ ' 



11 



is equivalent to CuciQ) on ly if the hemisphere- mixed sample is perfectly mod- 
elled, which implies N%$£(Q) /N^(Q) = 1. Figure || shows this ratio; the general 
agreement is good, although at Q < 0.8 GeV there is some indication that the 
ratio is below unity. 

Any important mismodeling of the Monte Carlo would indicate that the re- 
sults obtained with Cmc(Q) an d with C mix (Q) (since the Monte Carlo is used for 
normalising this correlation function) were not reliable. Therefore, the distribu- 
tion of Monte Carlo tracks in a sample free from BE correlations was compared 
to the same distribution in the data. Figure |6] shows the momentum spectrum 
of kaon candidate tracks of the hemisphere-mixed sample in the data and in the 
Monte Carlo normalised to the same total number of pairs. At low momentum, 
the Monte Carlo does not describe the data spectrum well and differences are seen 
at the ±15% level. Studies of the differential cross-section of kaons in hadronic 
events in both the data and the Monte Carlo events, generated with JETSET 
7.4 and tuned according to OPAL data, showed that the simulation predicted a 
kaon spectrum which is too soft. 

As a check of the stability of the results obtained with C m i X (Q), the Q spectra 
of both like-sign and hemisphere-mixed pairs were reweighted. Each pair of kaons 
in the Monte Carlo was reweighted by the product of the weights of each kaon 
in the pair, where the weight was obtained in bins of momentum by dividing the 
data momentum spectrum by the Monte Carlo spectrum. The final measurement 
of the correlation function did not change significantly after reweighting, A varied 
by +0.03 and Rq by —0.01 fm. The same exercise was done to check the results 
obtained with Cmc{Q)- 111 this case, A varied more significantly, by +0.12, and 
Rq by —0.01 fm. Thus, by the use of a double-ratio technique, the correlation 
function C mix (Q) was found to be less sensitive to the Monte Carlo mismodeling 
than the correlation function Cmc(Q)- 

5.6 Residual BE correlations in the reference sample 

As suggested in [fT5f1 , residual BE correlations could be a source of imperfection 
in the reference sample. To check that the hemisphere-mixed sample was free 
of effects induced by the BE correlations, a Monte Carlo study was done using 
the generator JETSET 7.4. Hadronic events were generated with and without 
BE correlations, with the BE correlations simulated assuming the Goldhaber 
parametrisation described in Section |[ The shape of the Q distribution of the 
hemisphere-mixed pairs remained unchanged after including the BE correlations 
in the generation. The ratio N™f x h BE \Q) / 7V™* W BE {Q) was consistent with 
unity in the full range of Q, demonstrating that the reference sample was free 
from effects due to residual BE correlations. 
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5.7 Sources of kaons 



In a substantial fraction of the kaon pairs, one or more of the kaons result from 
a long-lived particle decay — in such cases the kaon pairs cannot exhibit BE 
correlations. It is therefore useful to separate the various sources of kaons and 
to classify the parent particles to estimate the maximum possible value of A, 
as suggested in On the other hand, some studies p7^p9| have suggested 
that the correlation function is narrowed by the contribution of decay products 
of long-lived sources, and also that resonance decays induce a pair-momentum 
dependence of the radius. 

The kaon sources as predicted by the JETSET Monte Carlo simulation are 
given in Table [| These have been classified as in P,|27| into two main groups: 
long-lived sources with life-time cr > 10 fm, and short-lived sources with life-time 
cr < 10 fm. The table shows that the fraction of kaon pairs at low Q (< 0.6 GeV) 
with at least one kaon from a short-lived source is about 81%, so that the fraction 
of pairs in which both kaons arise from short-lived sources is about 66%. 

The pairs from short-lived sources cannot be identified in the data, and so 
the final results of this analysis were not corrected for the effect of short-lived 
sources because such a correction would be based on a Monte Carlo model with 
its inherent uncertainties. However, to illustrate the magnitude of the effect, a 
correction was applied to the correlation function using the estimated fraction 
of short-lived sources (66%), assuming that kaons from long-lived sources do not 
contribute to the correlations. The fitted parameters of the corrected correlation 
function are: A = 1.27 ± 0.31 and Rq = 0.55 ± 0.07 fm, where the errors are 
statistical only. 





origin of kaons 


fraction of kaons in 
pairs with Q < 0.6 GeV 


long-lived sources 


b, c hadron decays 


12% 


life-time cr > 10 fm 


0(1020) 


7% 


short-lived sources 


string fragmentation 


40% 


life-time cr < 10 fm 


fo(980) 


1% 




K*(892) 


27% 




other sources 


13% 



Table 3: Origins of kaons in like-sign pairs with low Q in the Monte Carlo. 



5.8 Final-state interactions 

Charged kaons are subject to both the Coulomb and the strong interactions. In 
principle, every pair of like-sign kaons from short-lived sources should be corrected 
for these interactions in the data (but not in the Monte Carlo, where they were 
not simulated) . To apply a correction to all pairs would result in an overestimate 
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of the value of the strength parameter [[30]. As in Section |5.7| , the final results of 
this analysis were not corrected for the electromagnetic and strong interactions 
because of the model-dependence of such corrections. 

As a check of the possible magnitude of any correction, the electromagnetic 
repulsion of like-sign pairs was corrected for in the data. The like-sign kaon pair 
Q spectrum was corrected using the Gamow factor |3]J G(r]) = 27irj j (e 27rr? — 1), 
where r] = a em m K /Q, a em is the electromagnetic coupling constant and m K the 
kaon mass. On the assumption that all pairs are from short-lived sources, the 
fitted strength and radius of the Coulomb corrected correlation function are: A = 
0.92 ± 0.25 and Ro = 0.61 ± 0.17 fm, where the errors are statistical only. The 
correlation function was also corrected for both the effect of short-lived sources 
as in Section ^T7| and the Coulomb effect, resulting in the fitted parameters: A = 
1.36 ± 0.55 and Rq = 0.58 ± 0.11 fm, where the errors are statistical only. 



6 Discussion and conclusions 

Bose-Einstein correlations have been measured in identified pairs of charged kaons 
in hadronic Z° decays using the OPAL experiment at LEP. The analysis was per- 
formed in events with a clear two-jet topology, a requirement which was neces- 
sary to obtain a suitable reference sample using a hemisphere-mixing technique. 
Monte Carlo simulation was used to correct for imperfections in the reference 
sample by use of a double-ratio for the correlation function. The enhancement 
was parametrised using a Gaussian formula, resulting in a strength 

A = 0.82 ± 0.22 t ° il 

and a kaon emitter radius 

R = 0.56 ± 0.08 1 g;gg fm. 

A definite conclusion from the present analysis is a confirmation of the re- 
sults of reference ]TTJ, that there are indeed BE correlations in K ± K ± pairs from 



hadronic Z° decays. This implies that there should be such correlations in K°Kg 
pairs [0 ; therefore it is unlikely that the previously observed threshold enhance- 
ments [jS^ITTH can be attributed entirely to the fo(980) decay into kaons. 

Values of A and Rq, as measured in hadronic Z° decays at LEP for various 



particle types, are listed for comparison in Table £|. Since there is evidence [g6 
that A and R may depend on the event topology, the table gives the type of 
event used in the various measurements. The reference sample types used in 
each analysis are also listed: these may be event- or hemisphere-mixed, unlike- 
sign or Monte Carlo pairs. 

In all events and for correlations measured using the unlike-sign reference 
sample, the radius of charged pion emitters varies between 0.8 and 1.0 fm while 
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pair 



A 

(stat. 



[sys.) 

U7T7 

0.12 

0.25 



Ro (fm) 
(stat.) 



0.06 



ref. 

sample 



events 



experiment 



K ± K 



0.82 ± 
0.82 ± 



0.22 I 
0.11 ± 



0.56 ± 0.08 



0.48 ± 0.04 ± 0.07 



mixed f 
unlike 



two-jet 
all 



this analysis 
DELPHI NTH 



7T ± 7T ± 



1.14 ± 
0.96 ± 
0.61 ± 



0.23 ± 
0.21 ± 
0.16 ± 



0.32 
0.40 
0.16 



0.76 ± 0.10 ± 0.11 
0.65 ± 0.07 ± 0.15 
0.55 ± 0.08 ± 0.12 



MC 
MC 
MC 



all 
all 
all 



OPAL 
ALEPH [| 
DELPHI Wl 



0.67 ± 
0.65 ± 
0.40 ± 
0.62 ± 
0.35 ± 
0.45 ± 



0.01 ± 

0.02 

0.02 

0.04 

0.04 

0.02 



0.02 



0.96 ± 0.01 ± 0.02 
0.91 ± 0.01 
0.49 ± 0.02 
0.81 ± 0.04 
0.42 ± 0.04 
0.82 ± 0.03 



unlike 
unlike 
mixed f 
unlike f 
mixed f 
unlike f 



all 

two- jet 
two- jet 
two- jet 
two- jet 
all 



OPAL P| 



OPAL |6| 
ALEPH [0 



ALEPH @ 
DELPHI ]TB| 
DELPHI H 



Table 4: Summary of the parameters A and Rq measured at LEP using the 
Gaussian parametrisation and for different types of identical pairs. The results 
marked with f were obtained using the double-ratio technique for the correlation 
function. All the it it results were corrected for Coulomb interactions and the 
ones marked with $ were in addition corrected for non-pion impurities. The 
uncertainties are statistical and systematic as labelled. 



the radius of the charged kaon emitters is 0.48 ± 0.08 fm. This gives the relation 
i?o( 7r±7r± ) > Ro(K ± K ± ) — a mass dependence of the emitting source, as already 
pointed out in 0, I n two-jet events, the measured radius of charged pion 
emitters is seen to vary between 0.4 and 0.9 fm, inconsistent within the quoted 
errors. This large variation is usually attributed to the choice of the reference 
sample: in the case of the unlike-sign reference sample the radius is about 0.8- 
0.9 fm; in the case of the mixed reference sample the radius is about 0.4-0.5 fm. 
The comparison of results obtained using the event- or hemisphere-mixing tech- 
niques and a double-ratio for the correlation function shows that the present 
measurement of the radius of the kaon source, 0.56 ± 0.11 fm, is compatible with 
the previous measurements of the radius of pion sources and does not support 
a strong mass dependence of the emitting source. However, both measurements 
of Rq for kaon pairs are considerably larger than that obtained for AA pairs, 



0.141 fm 
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Figure 1: Q distributions of like-sign and hemisphere- mixed KK pairs in 
the data (a) and in the Monte Carlo without simulation of BE correlations 
(b). The distributions were corrected for non-KK impurities. The hemisphere- 
mixed distributions were normalised to the like-sign distributions in the region 
0.6 < Q < 2.0 GeV, indicated by the arrows, where no BE correlations are ex- 
pected. The errors are statistical only. 
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Figure 2: Ratio N ++ (Q)/N mix (Q) for KK pairs in the data and the Monte 
Carlo. The first bin has double width due to the small statistics. The errors are 
statistical only. 
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Figure 3: C m i X {Q) BE correlation function in K ± K ± pairs and the fit using Equa- 
tion |] superimposed. The first bin has double width due to the small statistics. 
The errors are statistical only. 
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Figure 4: Cmc{Q) BE correlation function in K ± K ± pairs using Monte Carlo 
like-sign pairs as the reference sample and the fit using Equation f| superimposed. 
The first bin has double width due to the small statistics. The errors are statistical 
only. 
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Figure 5: Comparison of the hemisphere-mixed Q distributions in the data and 
the Monte Carlo using the ratio N^(Q)/N™£(Q). The errors are statistical 
only. 
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Figure 6: Momentum distribution of kaon candidates in the data and the Monte 
Carlo. The spectra for true kaons, pions and protons in the Monte Carlo are also 
shown. 
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